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Abstract 

A many body Hamiltonian involving the mean field for a projected spherical single particle 
basis, the pairing interactions for alike nucleons and the dipole-dipole proton-neutron interactions 
in the particle- hole (ph) and particle-particle {pp) channels is treated by the projected gauge 
fully renormalized proton-neutron quasiparticle random phase approximation (PGFRpnQRPA) 
approach. The resulting wave functions and energies for the mother and the daughter nuclei are 
used to calculate the 2vf3f3 decay rate and the process half life. For illustration, the formalism is 
applied for the decay ^"^"^Mo— )• ^"^"^Ru. The results are in good agreement with the corresponding 
experimental data. The Ikeda sum rule {ISR) is obeyed. The gauge projection makes the pp 
interaction inefficient. 

PACS numbers: 23.40.Hc,23.40.-s,21.10.Tg,21.60.Jz, 13.10,+q 
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Double beta decay is one of the most exciting topic of nuclear physics since the rate 
of the process is obtained by combining formalisms of electroweak interaction with those 
yielding nuclear matrix elements. Due to this feature it represents a sensitive test for both 
collaborating fields. The 2i//3/3 process is interesting by its own but is also very attractive 
because it constitutes a test for the nuclear matrix elements (m.e.) which are used for the 
process of Oi//3/3 decay. Discovery of this process may provide an answer for the fundamental 
question, whether neutrino is a Mayorana or a Dirac particle. The subject development is 
reflected be several review papers [l|-l7|. Our contribution described in this letter concerns 
the 21/ process, which can be viewed as two consecutive and virtual single /3~ decays. The 
formalism yielding closest results to the experimental data is the proton-neutron random 
phase approximation (pnQRPA) which includes the particle-hole {ph) and particle-particle 
{pp) as independent two body interactions. The second leg of the 2z//3/3 process is very 
sensitive to changing the relative strength of the later interaction, denoted hereafter by Qpp. 
It is worth mentioning that the two body interaction of ph type is repulsive while that of pp 
nature is attractive. Due to this feature there is a critical value for Qpp for which the first root 
of the pnQRPA equation vanishes. Actually, this is the signal that the pnQRPA approach 
is no longer valid. Moreover, the Qpp value which corresponds to a transition amplitude 
which agrees with the corresponding experimental data is closed to the mentioned critical 
value. That means that the result is not stable to adding corrections to the RPA picture. 
The first improvement for the pnQRPA was achieved by one of us (AAR) in Refs. 8|, l9|, by 
using a boson expansion (BE) procedure. Later on another procedure showed up, which 
renormalized the dipole two quasiparticle operators by replacing the scalar components of 
their commutators by their average values Such a renormalization is inconsistently 

achieved since the scattering operators are not renormalized. This lack of consistency was 
removed in Ref. 
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12l | where a fully renormalized pnQRPA is proposed. 
Unfortunately, all higher pnQRPA procedures mentioned above have a common drawback 
of violating the Ikeda sum rule (ISR) by an amount of about 20-30% [13]. It is believed 
that such a violation is caused by the gauge symmetry breaking. Consequently, a method 



of restoring this symmetry was formulated by the present authors in Ref. 

In this paper the results of Ref. [l^ are improved in two respects: a) aiming at providing a 
unitary description of the process for the situations when the nuclei involved are _spherical or 
deformed, here we use the projected spherical single particle basis defined in Ref.[l5| and used 
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for double beta decay in Ref . [16|, ll7|] . b) the space of proton- neutron dipole configurations 
is split in three subspaces, one being associated to the single (3~, one to the (3~^ process, and 
one spanned by the unphysical states. A compact expression for the dispersion equation of 
energies is obtained from the linearized equation of motion of the basic transition operators 
corresponding to the two coupled processes. The numerical application is made for the 2i>/3/3 
process ^°°Mo — j.^'^o Ru. Aiming at a selfcontent presentation, minimal details are necessary. 
According to Ref. 15| the projected spherical basis is defined as: 

Kil (d) = KiAnllr^^jI)^,] ^ Ki.Kil id), (1) 

where Pmk denotes the angular momentum projection operator, \nljm) is the spherical shell 
model state and is an axially deformed coherent state describing the ground state of a 

^g = exp[d{bt,-b2o)]\0)b. (2) 



phenomenological core in terms of quadrupole bosons fcL? ^: 



Here |0)f, denotes the vacuum state for the quadrupole bosons. The set of functions (d) is 
an orthogonal single particle basis. The single particle energies e^;^- are obtained by averaging 
a particle-core Hamiltonian with the corresponding basis states. In order to keep close to 
the Nilsson model, where on each Q state one can distribute 2 nucleons, here we change the 
norm of the projected states such that this restriction holds also for each / state: 

($^*'i<^>n = i^E("^a"'i<fn=2. (3) 

M 

Thus, the wave functions used to calculate the m.e. should be multiplied with the statistical 



factor ^2/(2/ + 1). 

We suppose that the states describing the nuclei involved in a 2z//3/3 process are described 
by a many body Hamiltonian which may be written in the projected spherical basis as: 

2 t ^ t 

H = E 21+1 ~ ^ra)CrcxIM^raIM " X! '^PralP'ral' 

T,a,I,M T,a,I,l' 

+ 2x E (3;ipn)(3l,{p'n'){-r -2x1 E P-{pn)P^,{p'n'){-r, (4) 
/ / / / 

pn]p n pn;p n ;/x 

where c\^jj^^{craiM) denotes the creation (annihilation) operator of one nucleon of the type 
r(= p, n) in the state with a being an abbreviation for the set of quantum numbers nlj. 
The Hamiltonian H contains the mean field term, the pairing interaction for alike nucleons 



and the Gamow- Teller dipole-dipole interaction in the ph and pp channels, characterized by 
the strengths x and Xi? respectively. Passing to the quasiparticle representation through 
the Bogoliubov-Valatin transformation: 



\I-M 



r = p,n, U^j + V,^, = l, 



(5) 



the first two terms of H are replaced by the independent quasiparticles term, J2 E 

while the ph and pp interactions are expressed in terms of the dipole two qp and the qp density 

operators: 



(6) 



In Ref. Ill], we showed that all these operators can be renormalized as suggested by the 
commutation equations: 



Bl{k),Al,{k')] ^ [Bl{k),A^,>{k') 



0, 



p p 

n p 



k 



(7) 



Indeed, denoting by CfJ^^ and Cf^j^ the averages of the right hand sides of d?]) with 
the renormalized RPA vacuum state, the renormalized operators defined as Ai^{k) = 
^==^41^, Bi^{k) = —j==Bif^, obey boson like commutation relations: 



ic 



(1) 



c 



Ai^{k),Al ,{k')] = 6k,k'S, 



'If,- 



Bi^{k),Bl^,{k')\ = 6k,k'Sf,,f,'fk, fk = sign{Cf^). 
Further, these operators are used to define the phonon operator: 
Cl = Y: [x{k)Al{k) + Z{k)Dl{k) - Y{k)A,.,m-f-^ - W{k)D,.,{k){-Y-^ 



(8) 



(9) 



where D\^{k) is equal to B\^,{k') or Bi^{k) depending on whether fk is + or -. The phonon 
amplitudes are determined by the equations supplied by the operator equations: 



(10) 



Interesting properties for these equations and their solutions are discussed in our previous 



pubhcations 
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12| . Here we mention one of these features. The renormahzed ground state 



is a superposition of components describing the neighboring nuclei (A^ — 1, Z + 1), (A^+1, Z — 
1), (A^ + IZ + 1), {N — 1, Z — 1). The first two components conserve the total number of 
nucleons (N+Z) but violates the third component of isospin, T3. By contrast, the last two 
components violates the total number of nucleons but preserve T3. Actually, the last two 
components contribute to the violation of the ISR. One can construct linear combinations 
of the basic operators A"^ , A, B"^ , B which excite the nucleus (A^, Z) to the nuclei (A^ — 1,Z + 
1), (A^ + 1, Z — 1), (A^ + 1, Z + 1), (A^ — 1, Z — 1), respectively. These operators are: 

•4^.{pn) = UpVnAl^ipn) + UnVpA,^.^{pn){-)'-^' + UpU^BlXpn) - K.K5l,_^(pr^)(-)l-^ 
Ai^ipn) = U^V^A.^ipn) + UnVpAl_^{pn){-y-^ + UpUnB^,,{pn) - 

Al^M = U,UnA{^{pn) - V,VnA^^^^{pn){-Y-^ - U,VnB\^{pn) - V,f/„fii,_^(pn)(-)i-^ 
Ai^M = U,UnA^,{pn) - VjyrAi-M)^-^''' - UpV^B.^ipn) - VpUnBl^^ipn)^-)'-^ . 
Indeed, in the particle representation these operators have the expressions: 



^p^n 



^p^n 



In terms of the new operators the many body Hamiltonian is: 



H 



5Z EjrjO'rjm'^Tjm 



^pn;p'n' 



(^pn;p'n' 



^1 11' 



{h\\^\\^n){Ip'\W\\In')- 



Here Eri denotes the quasiparticle energy. Since we are interested in describing the harmonic 
modes which preserve the total number of nucleons, we ignore the Xi term. Indeed, this 
term defines a deuteron type excitation, which will be studied in a separate publication. 
The equations of motion of the operators defining the phonon operator are determined by 
the commutation relations: 



Ai^{pn),A\^,{p'n 



Jn 



f/2 - -1/2 ^ 



Jp 



(12) 



The average of the r.h. side of this equation with the PGFRpnQRPA vacuum state is 
denoted by: 

1 .... .... , 1 



Diipn) = U' 



^ n ' 



2/„ + l 



2h + l 



iU'-V'){N,). 



(13) 



(11) 



The equations of motion show that the two qp energies are renormahzed too: 

E^^^ipn) = E,{Ul - V^) + E^{V^ - Ul). (14) 

The space of pn dipole states, <S, is written as a sum of three subspaces defined as: 

5+ = {(p, 71)1^1 (pn) >0, >0,}, 5_ = {(p, n)|L>i(pn) < 0, E'^'^ipn) < Q,} , 

= <S-(5+ + 5_). (15) 

In 5+ one defines the renormahzed operators: 

- 1 - 1 

At^ipn) = -^==A\i,{pn), Ai^,{pn) = Ai^,{pn), (16) 

yDiipn) yDiipn) 

while in <S_ the renormahzed operators are: 

•^iVW = -J===Al^,ipn), J'i^(pn) = ^ ^ A\^ipn). (17) 
Vl^iWI \/\Di{pn)\ 

Indeed, the operator pairs Ai^, A\i^ and J^i^, satisfy commutation relations of boson 
type. An RPA treatment within Ssp would yield either vanishing or negative energies. The 
corresponding states are therefore spurious. 

The equations of motion for the renormahzed operators read: 

H,Al] ^ E^^^{pn)Al + 2x E 

(pini)e5+ 

(pini)G5_ 

'H,n,] = \E^^^pn)\n^ + 2x E ^inU-^l 

(pini)e5_ 

+ 2X E ^linnAi-,{-l)'-', 

(pini)e5+ 

H,A 



5 -^IfJ. 



where: 



: -E^'^-{pn)A,, - 2x E '^Sp.n.^i. 

(pini)eS+ 

(pi?ii)e5- 

H,f^,] = -\K^-{pn)\T,,-2x E ^Ui-^iM 

(pini)G5_ 

- 2X E ^Sp.n.4,-,(-l)'-^ (18) 

(pini)g5+ 



The phonon operator is defined as: 



(20) 



with the amphtudes determined by the equations: 



1 ifj, 



(21) 



The compatibihty condition for the homogeneous system of equations determining the 
phonon amphtudes yields two dispersion equations for u: 



2x [nt - n-) = 1, 2x (-7^+ + 7^:) = i, 



with 



7^1 



E 



T, 



p\n\ 



(pini)e<SH 



UJ ± E'"^"(pini; 



-, 7^: 



E 



(pini)e5- 



(22) 



(23) 



The phonon amphtudes can be analytically determined: 



X{pn) = 2x 



- pn 



C, W{pn) = -2x 



rp2 

pn 



Z{pn) = 2x- 



T? 



pn 



-C\ Y{pn) = -2x 



pn 



'UJ - |E^^"(pn)| 
The constant factors C and C' have the expressions: 



UJ + E-''''^{pn) 



(24) 



= Ax' 



nt-n- 



C'~' = 4x' 



TT -n] 



(25) 



In order to solve Eqs.f l22p we need to know Di{pn) and, therefore, the averages of the gp's 
number operators, Np and Nn- These are written first in particle representation and then 
the particle number conserving term is expressed as a linear combination of A'^A and J-'^J-' 
chosen such that their commutators with A^ , A and J-"^, J-" are preserved. The final result is: 

{Np) = V;{2Ip + l) + 3iU'p-Vp')i Y: D,ip,n)iY,{p,n)Y- ^ D^ip,n){W,{p,n)r), 

n ,k n ,k 

{Nn) = V^{2I,, + 1) + 3{U'^-V^){ Y: D,{p,n){Y,{p\n)r- Y D,ip\n){W,{p ,n)Y). 



P 5^ 

(p',n)eS+ 



p 1^ 

{p',n)eS- 



Eqs. ( l22l) . (1261) are to be solved iteratively. It is worth mentioning that using the quasiparticle 
representation for the basic operators J-"/^, ^i^_^(— 1)^^'*, J-'i^_^(— )^~'^, one obtains for 
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an expression which involves the scattering pn operators. Thus, the present approach is 
indeed the PGFRpnQRPA. 

The formahsm presented above was used to describe the 2z//3/3 process. If the energy 
carried by leptons in the intermediate state is approximated by the sum of the rest energy 
of the emitted electron and half the Q- value of the double beta decay process 

AE =^Qpp + mec\ (26) 

the reciprocal value of the 2z//3/3 half life can be factorized as: 

(T,7f)-i = F|McT(0+-^0+)p, (27) 

where F is an integral on the phase space, independent of the nuclear structure, while Mgt 
stands for the Gamow- Teller transition amplitude and has the expression : 

In the above equation, the denominator consists of three terms: a) AE, which was already 
defined, b) the average value of the k-th PDFRpnQRPA energy normalized to the particular 
value corresponding to k=l, and c) the experimental energy for the lowest 1"^ state. The 
indices carried by the f3~^ operators indicate that they act in the space spanned by the 
PGFRpnQRPA states associated to the initial (z) or final (/) nucleus. The overlap m.e. 
of the single phonon states in the initial and final nuclei respectively, are calculated within 



GPFRpnQRPA. In Eq. (!28|) . the Rose convention for the reduced m.e. is used [18 |. 

Note that if we restrict the pn space to 5+, Mgt vanishes due to the second leg 
of the transition. Indeed, the m.e. associated to the daughter nucleus is of the type 
/(0|(cJjCp)i^(cJjCp)i^|o)/, which is equal to zero due to the Pauli principle restriction. Also, 
we remark that the operator plays the role of a /3~ transition operator, while when 
is applied on the ground state of the daughter nucleus, it induces a f3~^ transition. Therefore, 
the 2/3 decay cannot be described by considering the /3~ transition alone. 

For illustration, we present the results for the transition ^°'^Mo— )'^°'^Ru. For this case the 
energy corrections involved in Eq.( l28l) are: 



AE = 2.026MeV, Ei+ = O.OMeV. (29) 
The parameters defining the single particle energies are those of the spherical shell model, 




FIG. 1: The single /3~ strength for the mother nucleus, ^'^'^Mo (left panel), and the strength for 
the daughter nucleus, ^^''Ru, folded by a Gaussian with a width of O.SMeV, are plotted as functions 
of the corresponding energies yielded by the present formalism. 





d k Gp[MeV] G„[MeV] ISR log ft x[MeV] 


iooMo 


-1.5 6.1 0.18 0.263 16.01 ioomo^^/^^^ooTc 0.207 

4.5 5.29 




-0.6 6.1 0.18 0.25 12.01 i°°Tc^i°°Ru 0.207 

4.66 5.41 



TABLE I: The deformation parameter d, the pairing interaction strengths for protons (Gp) and 
and neutrons (G„) and the GT dipole interaction x used in our calculations. We also give the 
parameter k relating the quadrupole coordinates and bosons (this is involved in the expression of 
the single particle energies) as well as the resulting log/t values characterizing the j3~^ /EC and /3_ 
transitions of ^"^Tc. The results for log/t values, given in the right column, are compared to the 
experimental data from the left column. 



the deformation parameter d and the parameter k relating the quadrupole coordinate with 
the quadrupole bosons: 

«2;. = 4(^1-(-1)"Vm)- (30) 



These are fixed as described in Ref. 



17|. The core system is defined by {Z,N) = (20,20). 



Labeling the states according to their energies ordering, the single particle space is defined 
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by the indices interval [11,62]. The dimensions for the spaces (iS+,iS_) are (163,1) and 
(158, 5) for the mother and daughter nuclei, respectively. The dimension for S is 208 for 
both the mother and the daughter nuclei. The strength of the dipole pn was taken to be 

^ = l^^e^- (31) 
This expression was obtained by fitting the positions of the GT resonances in ^°Ca, ^"Zr and 



208 



Pb [19|]. The value obtained for the Mo and Ru isotopes is that given in Table 1. 
Using these input data we calculated the distribution of the strengths with the re- 
sult shown in Fig.l. The energy intervals where both distributions are large, contribute 
significantly to the double beta transition amplitude. 

Calculating first the GT transition amplitude and then the Fermi integral with Ga = 
1.254, as in Ref.Q], we obtained the result: 

\Mgt\ = 0.234, Ti/2 = 7.89 ■ lO^V- (32) 



This result should be compared with the experimental results 20|, l21| : 



Ti/2 = (8.0 ± 0.16) ■ 10 V, Ti/2 = {O.Ubt^om) " lO'V (33) 

Another experimental result concerns the summed strength for the /3~ transition: Bgt- = 
26.69. Quenching the theoretical result by a factor 0.6, as to account for the missing strength, 
one obtains the value of 28.9. The intermediate odd-odd nucleus, ^"^^Tc, can perform the 
transition P~^/EC, feeding ^"^^Mo, or the /3~ transition to ^°'^Ru. The measured log ft values 
for these transitions, are given in Table 1. The corresponding theoretical results are obtained 
by means of the expression: 

fe = !!!° (34) 



In order to take account of the effect of distant states responsible for the "missing strength" 
in the giant GT resonance [4] we chose qa = 1-0. 

Summarizing the results of this paper, one may say that restoring the gauge symmetry 
from the fully renormalized pnQRPA, one obtains a realistic description of the transition 
rate and moreover the ISR is obeyed. As shown in this paper, it seems that there is no 
need to include the pp interaction in the many body treatment of the process. 
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